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Abstract 

Decoherence effects associated to the damping of a tunneling two-level system are 
shown to dominate the tunneling probability at short times in strong coupling regimes 
in the context of a soluble model. A general decomposition of tunneling rates in dissi- 
pative and unitary parts is implemented. Master equation treatments fail to describe 
the model system correctly when more than a single relaxation time is involved. 

1 Introduction 

Motion in tunneling phenomena can be seen quite generally as resulting from dispersive 
effects in the time evolution of spatially localized quantum states which are not energy 
eigenstates. Since, from this point of view, quantum interference plays a central role in these 
phenomena, one may ask how are they affected by quantum decoherence processes which 
occur as a result of coupling to other degrees of freedom. Effects of additional degrees of 
freedom in tunneling processes have been studied now for a long time and in a variety of 
contexts []T] , but emphasis falls usually on the question of the observable changes on inclusive 
properties such as tunneling probabilities. However, these questions are often addressed 
in semi-classical terms, including generalizations of the WKB approach [0, which tend to 
becloud the dispersive character of the quantum dynamics. 

In order to analyze the dispersive dynamics of tunneling we consider in this paper the 
simple but not unrealistic case of two-level approximations to symmetric, bound, bi-stable 
systems whose energy spectrum is characterized by a structure of doublets. Each member 
of one doublet consists of symmetric/antisymmetric superpositions of localized states which 
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may be seen as coupled by effective tunneling amplitudes. In the two-level approximation 
this is conveniently expressed in terms of a spin- 1/2 algebra. Dissipation and decoherence 
phenomena can be introduced in this description by coupling the two-level system to another 
system, possibly with a continuous energy spectrum. A collection of harmonic oscillators is 
frequently used in this connection, leading to the "spin-boson" model analyzed by Leggett 
and collaborators |J. There the effects of dissipation on the tunneling rates are studied 
independently of explicit control of properties of the tunneling state relating to decoherence 
processes. A simpler, soluble model, which allows for a control of this sort, results when 
the two-level system is suitably coupled to a generic additional system with continuous 
spectrum, as described in detail in section |3| below. In this model the coupled dynamics 
can be reduced to the problem of the spreading of a product "doorway state" due to its 
coupling to a continuous background also of product states, which allows for a treatment in 
terms of the techniques developed long ago by Fano [|J in the context of atomic scattering 
theory. While this schematic model is less "realistic" from the point of view of energy-loss 
mechanisms, it retains the overall dynamical features of the spin-boson model relating to 
coherence-loss, which is our main concern here. The model allows for the calculation of the 
reduced density associated with the two-level system at all times. Model dependences are 
indicated by differences in the results obtained in the two cases. We find, in particular, that 
transition rates out of a localized initial state can be substantially enhanced by decoherence 
in strong coupling regimes. 

The use of Master equations for the description of damped two-level systems is well 
known || . The usual derivation of these equations || assume weak-coupling to an harmonic 
oscillator reservoir and high temperature, so that their use in strong coupling and zero tem- 
perature situations is not warranted. Comparison of the exact model results with solutions 
of these Master equations agree for particular forms of the coupling independently of the 
overall coupling strength and at zero temperature. In general, however, the perturbative 
ingredients of the Master equation derivations appear as limiting factors. 



2 Dissipative and unitary tunneling rates. 

The quantum state of a two-level subsystem can be generally described in terms of a (time 
dependent) hermitian reduced density operator of unit trace p(t). The simple two-level 
effective Hamiltonian 

Ha = | (CT 3 + 1 CT ), 

where l a stands for the unit 2x2 matrix, can be understood as describing the tunneling 
between the "localized" ("left" and "right") states 

10 - ' +) + '" } and |r> ee l+) ~ H 

where |±) are the eigenvectors of the diagonal Pauli matrix 03, and hence also energy eigen- 
states. The time dependent density matrix corresponding to the initial condition 
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p(t = 0) = \l)(l\ (1) 
is given in the |±) basis as (with % = 1) 



p(t) = U e L el ^=p 2 (t). (2) 

As is well known, the property of idempotency is equivalent to the purity of the quantum 
state described by pit). A tunneling probability P(t) can be defined as the probability 
of finding this system in the complementary localized state |r) at time t, and is readily 
calculated as 

P(t) = Tr[|r)(r|p(t)] = sin 2 |. (3) 

When the two-level subsystem is coupled to additional degrees of freedom its state at time 
t is still described by a reduced density matrix which is non-negative, hermitian and of unit 
trace, but in general not idempotent. In fact, the quantity 5{t) = 1 — Tr[p(t) 2 ] (referred to as 
the linear entropy or the idempotency defect) is often used as a measure of the decoherence 
suffered by the two level subsystem through its coupling to the other degrees of freedom. 

A most convenient way of representing such a reduced density matrix makes use of its 
eigenvalues and eigenvectors 

p(t)\t,k) = Pk (t)\t,k), fc = l,2; (t,k\t,k') =8 W , 0< Pk (t)<l 
in terms of which one can write 

p(t) = Y,\t,k)p k (t)(t,k\. (4) 

k 

It should be noted that in general both the eigenvectors and the eigenvalues in (Q) are 
time-dependent. While the eigenvectors evolve unitarily (as they define a time-dependent or- 
thonormal basis in the quantum phase-space of the two-level subsystem) the time-dependence 
of the eigenvalues p k (t) reveal the change in time of the coherence properties of the state of 
the two-level subsystem. The idempotency defect is in fact given in terms of the eigenvalues 
Pk(t) as 

5=l-J2Pk(t) 2 = 2p 1 (t)(l-p 1 (t)) (5) 
k 

where use has been made of the unit trace property of p(t). In the case of eq. (Q), it can be 
easily checked that p± = 1 and p 2 = at all times, so that 8 = 0. Assuming that the reduced 
density is that which evolves from the initial condition ([!]), the tunneling probability is now 
given by 

P(t)=J2p k (t)\(r\t,k)\ 2 . 

k 
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In view of what has been said concerning the time-dependence of the ingredients of this 
expression, it is immediately clear that the tunneling rate 

k k ai 

splits into two contributions which can be unambiguously ascribed to changes in the coher- 
ence properties of the state of the subsystem and to its unitary evolution, respectively: 

Rd(t) = E Pk(t) I (r |t, k}\ 2 = (2\(r\t, 1) | 2 - l) Pl (t) (6) 
k 

and 

Ru(t) = T^Pk(t)j t \(r\t,k)\ 2 = (2 Pl (t) - 1) ±\(r\t, 1)| 2 . (7) 

The last forms of R u and Rd make use of the unit trace of p and also of the normalization of 
|r). In the case of the purely unitary time evolution under the simple two-level Hamiltonian 
H a one has of course Rd(t) = and R u (t) = R(t) = (e/2)sin et. This decomposition of 
transition rates in dissipative and unitary parts can be simply extended to cases involving 
reduced densities of larger dimensionality, as indicated in Appendix A. 



3 Two-level tunneling with coupling to a continuum. 

A simple model allowing for decoherence effects in the two-level tunneling processes is one in 
which the two-level subsystem described by the Hamiltonian H a is coupled to an additional 
"nondescript" set of degrees of freedom (which will be identified by means of a subscript 
b) to which is associated an energy spectrum containing one discrete state represented by a 
normalized state vector |0&) in addition to a continuum of states normalized in energy 
as (77 1 7/) = 5(r) — T)'). The model is characterized further by the Hamiltonian 

H = ^(<T 3 + i*)®i b + (\0 b )e (0 b \+ J\r) ® 

+\ f V dv \g(v)\v)<r-(0b\+g*(T,)\0b)<T + (v\] (8) 

Z Jr)o 

+ \ f V dr) W(v)\v)<r+{0b\+!/*(ri)\0b)<T-{v\] 

Z Jr)o 

where g(r/) and g'(r)) are coupling matrix elements (— r}\g\ +0&) and (+ri\g'\ — Ob) respectively, 
possibly dependent on rj. In the special case g(r)) = g'(rj) the coupling terms reduce to 

v* f V dri [\rj)g(rj)(O b \ + \O b )g*(rj)(rj\} 
Jvo 
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which resembles in its structure the coupling term of the spin-boson model 0. 

A basis in the quantum phase-space of the composite system described by equation (§) 
can be constructed as {|s0&), {|sr/)}}, s — +,—, rjo < r\ < fj, and the adopted structure 
of the coupling terms shows that the subspaces spanned by the two sets of basis vectors 
{|+0b), {| —??)}} and {| —Ob), {| + v)}} are closed under the action of H. Different choices for 
the model parameters e, e , the span of the continuum band rj to fj and the ^-dependence of 
the coupling functions g and g' allow for a considerable variety of dynamical regimes. Closest 
to the situation described by the spin-boson model is that in which eo = 1] = with g(r]) = 
g*(rj) = g'(i]). Some aspects of the relationship of these two models are briefly discussed in 
Appendix B. The particular version of the model which will be used in the following sections 
differs from this, however. A manifest feature of this version is the assumption that one has 

r]o e and e + e fj 

so that the two product states | ± Ob), with eigen-energies eo and eo + e fall well within a 
sufficiently wide continuum band of product states | ±77). This feature will in fact allow for the 
continuum spreading of the states | + Ob) and | — 0&) through their coupling with states | — rj) 
and I + rj) respectively with spreading widths appreciably larger than the doublet splitting 
e, a situation which will be seen to characterize strong damping regimes. Furthermore, since 
these two spreadings are dynamically independent, as they take place in different invariant 
subspaces of H, there will be in general two distinct relaxation times associated with the 
damped two-level subsystem. 

General state vectors in each of the invariant subspaces of H can be written as 

\^ a ) = ao\ +0 b ) + [ dr]A(r])\ -77) 
\*k) = b Q \-0 b )+ [ V d V B(rj)\+ V ). 
In terms of these expansions, the eigenvalue problem for H 

(E-H)\^) = 

can be reduced to two independent pairs of coupled equations with similar structure for the 
expansion coefficients 

(E-e-e )a { E) = f d V g*( V )A^( V ), 

(E-r,)A^( V ) = g( V )aP (9) 

and 
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(E-e )bi E) = / d V g'% V )B( E \ V ), 
(E-e-r])B^(r]) = g'(v)bP. (10) 



The solutions of these coupled equations has been given long ago by Fano[Q and van 
Kampen|7]]. In the case of eqs. @ one has 

, a (*), 2 = \mi ai) 

1 ° ' ~ [E-e-e G -F{E)¥ + ^\g{E)\± 1 ) 

with F(E) given by the principal value integral 

F{E)=V r^J^£. (12) 
J no E -i] 



The continuum coefficients are 

1 9(v)a { o E) (13) 



with 



P + z{E)5{E - rj) 



z(E) 



E-r) 

E — e — e — F(E) 



\g(W 



( pi 

The value of the amplitude % involves a phasing convention which is subsequently carried 
over to the A^ E \rj) through equation (|T3|). The solution of eqs. ( fTPf ) is of course completely 
analogous and does not have to be given explicitly here. 

Strictly speaking, when the integration limits of the continuum variable r] are finite, the 
spectrum of H may include discrete states outside the continuum range, subject to the in- 
dependence of the coupling parameters g(r]). Whenever needed, the discrete eigenvalues of 
H can be obtained (e.g. in the case of eqs. (|])) as solutions Ed < r)o or E4 > fj of the 
dispersion equation 

E — e — e = / dr\ 



'/o 



E-rj 



In this case eq. ([H|) is replaced by 



E d -v 



and is determined from the normalization condition (up to an arbitrary phase factor) 

as 
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, , -1 



(IE j 



E=E d 



The low-energy discrete solutions are specially relevant for comparison of this model with 
other models of the damping mechanism, such as the many-oscillators model (see Appendix 
B). 



3.1 Time evolution of a localized initial condition. 



Using the stationary states of H as determined above one can write the time evolution of 
the localized initial state 



\t = 0) 



!+> + !-> 

y/2 



®|o 6 ) = |0®|o 6 ). 



(14) 



Note that in the absence of coupling, g(rj) = g'{rj) = 0, the time evolution reduces to that 
found in section |^. In general, the component states | + Ob) and | — 0&) can be written in 
terms of the stationary states determined from the solutions of eqs. (0) and (|10D, so that 



\t) = e~ im \t = 0) 



-iHt 



vo 



When the continuum range extends over a sufficiently broad energy interval on both sides of 
e the discrete state amplitudes a^ d ^ and become very small and can be neglected. In 
what follows we assume this to be the case in order to avoid unessential complications. By 
re-expressing the stationary states in terms of the factorized bases, forming the total density 
\t)(t\ and taking its trace over the 6-states one obtains for the reduced density of the two 
level system at time t 



P(t) 



P++(t) p+_(t) 
p*{t) l-p++(t) 



with 



P++{t) 



1 
2 



1 + 



f dE\af ] \ 2 e^ Et 


2 


f dE\b[ E) \ 2 e- iEt 


') 











(15) 



and 
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p+-{t) 



1 

2 
+ 



J rjo J r) 



i(E-E')t 



dr]B( E \r))A( E '>{r)) 



(16) 



The last term of eq. (|T6| ) is the only one in which the continuum amplitudes A^ E '\rj) and 
B^ E \ri) intervene explicitly. This term describes a re-correlation of components of the state 
of the two-level system in the single continuum subspace of system b, to which they couple 
through g(t]) and g'{rf). In fact, a simple variant of the model that avoids this re-correlation 
process can be devised by adding a second continuum {|C)} orthogonal to the first, i.e. 
(?7|C) = 0, to which the component | — Of,) is coupled through g'((). It is easy to see that the 
resulting expressions for p ++ and p + _ are in this case identical to eqs. (|l^) and (|16D except 
for the absence of the last term in the latter. 

The tunneling probability at time t of the two- level subsystem, from the initial state |Z) 
to the complementarily localized state |r), defined in terms of the reduced density matrix 
pit) as in eq. (|]), is generally given as 



P(t) = Tr[|r)(r|p(t)] = -(l-2Rep + _) 
and the idempotency defect (|) of the reduced density is 

S(t) = 2 \ P++ (t) (1 - p ++ {t)) - \p+-(t)\ 2 } = 2detp(t). 



(17) 



3.2 Special case: single relaxation time. 

An interesting particular case is that in which g'(f]) = 0, so that | — 0?,) is a stationary state 
of H while | + 0&) is spread through its coupling to the | — rf) continuum. In this case the 
reduced density matrix elements are simply given in terms of the spectral distribution |a[) E ^ 2 , 
eq. (|TT|), of the state | + 0&) as 



P++(t) 



dE\a ( E) \ 2 e- tEt 



'/ii 



and 



A Jr)o 



-i(E-e )t 



g'{rj) = 



g'iv) = o. 



(19) 



(20) 



These expressions can be immediately evaluated in closed form when the matrix elements 
g{rj) are independent of i], and the energy shift F(E), eq. (|l"2D is slowly energy dependent 
over an interval of the order of 7t\g\ 2 in the neighborhood of E = e + e^. This will in fact be 
the case when the continuum range extends over sufficiently broad intervals both above and 
below this energy. Equation (|TTD is then well approximated by the Breit-Wigner form 
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i (£),2 W = J_ r 

where e# — e — eo — F{e R ) = and T = 2tt | 1 2 is the usual width parameter of the Breit- 
Wigner distribution, given in Golden Rule form. Note that for large fj — e$ ~ eo — ?7o one 
has F(e R ) ~ ^(eo) — so that e# — eo — e. Extending the integration limits to ±oo one 
obtains by simple contour integrations 

p++ (t) = \e- Tt and p+_(t) = I e -^- e °)*-ir (21) 



Pit) = -(l- e~¥ cos{e R - e )t) . (22) 



2 

so that 



2 

It should be kept in mind that the model assumptions allow T > e, in which case the initial 
time-dependence of P(t) is dominated by the decaying exponential factor. The decoherence 
process undergone by the reduced density as measured by the idempotency defect (||) is 
given by 

5(t) = le"" (l - e-^) (23) 
so that the tunneling and the decoherence rates 



and 



P(t) = ( - cos(e R - e )t + (e R - e ) sin(ei? - e )t 



5(t) = T -e~^ (2e" r * - 1 



2 

are both proportional to the damping width T, related to the single relaxation time defined 
as Td = 1/r. Moreover, P(t) contains the period of unitary evolution t u = 27r/(e R — eo) — 
27r/e as an additional, independent time scale. In the strong coupling regime <^ t u eq. 
(f22|) shows dominance of coherence loss effects in P(t) for times t < Td leading eventually to 
saturation at P = 1/2. The sorting of dissipative and unitary contributions to the tunneling 
rate P(t) is achieved through the calculation of the rates Rdit) and R u {t) introduced in 
section || As shown in eqs. (||) and (0), this involves the time dependence of the eigenvalues 
and eigenvectors of the reduced density. The relevant ingredients can be expressed in a 

straightforward way in terms of the reduced density matrix elements p ++ (t) and p H (t). 

Since the resulting expressions are lengthy and not particularly illuminating they will not be 
given here. Numerical results illustrating these features, obtained using the matrix elements 
given in eq. fl2~ip, are shown in figs. 1 and 2 respectively for the weak and strong coupling 
regimes. The figures show P(t), 5(t) and the tunneling rates P(t), Rd(t) and R u {t). The 
correlation of 5(t) with Rd{t) and the dominance of the latter over R u {t) in the strong 
coupling regime are clearly seen. 
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3.3 General case: two relaxation times. 



When both g(r]) and g'(rj) are different from zero the matrix elements of the reduced densities 
are given by the complete expressions (|TJ) and (|T6|) . Using the tunneling initial condition 
of subsection |3J] and under the assumptions used in subsection |37| the first of these can be 
evaluated to yield 



P-\ 



.w - \ (i 



,-r't 



(24) 



where T', defined as the width parameter V = 27r|g'| 2 (g' having been assumed to be inde- 
pendent of rj) of the distribution 



\b {E) \ 2 



r 



2tt (E - e' R f + 



r' 2 



introduces an additional time scale in the evolution of the composite system. As for p + _ one 
gets, from eq. (|i~6|), using eq. ( |13D and the methods of ref. |4]] 



-i(e R -e' R )t- 



r+r' 2VTT 7 



+ 



r + p 



i(e R -e' R )t { 1 _ e -^t 



(25) 



where the last term is due to the re- correlation processes mentioned in subsection |3.1|. Note 



that under conditions in which the shift functions -F(e^) and F'(e' R ) are small one has 
e R — e R — e - This will be assumed to be the case from here on in order to save notation. 
The tunneling probability is given by eq. ([D]) as 



P(t) 



1 cos et 

2 2~ 



2 V / TT 7 

r + F 



+ e~ 



r+r' 



1 



2 v / TT 7 ' 

r + r 



(26) 



which shows that the re-correlation processes give rise to steady state oscillations of P(t) 
after the transient associated with the decaying exponential has died down. The amplitude 
of the steady-state oscillations depends on the damping widths T and T', however, and is 
maximum when T = V, in which case the transient actually cancels out and the undamped 
form of P(t), eq. (||), is recovered. 

This last result does not imply, however, that there is no decoherence when T = V. The 
general expression for 5(t) obtained from eq. (0) can be written as 



5(t) 



1 + e~ rt - e- r,t 



2Vrr 
~r + r 

and reduces, when T = T', to 



1 — e 2 



1 - e- rt + e- T,t 



r+r j 

e 2 cos 2et + 



-(T+r')t 



r + r 



r+r' 



1 — e 2 



5(t) 



r=r' 



2e- r Ml-e- ir )sin^ et 



rn -2 
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which vanishes only asymptotically, as t — > oo. In the general case r ^ T' with T,T' ^ 0, 
one has the asymptotic limit 



The restriction r, V 7^ is necessary since the limit t —>■ oo does not commute with the 
limit T' — > 0. If the alternate version of the model with two orthogonal continua is used, all 
terms involving factors VTT'/ (T + V) are absent, both in P(t) and in 5(t). In this case the 
two-level system decoheres maximally when t — > oo as two orthogonal components of the 
two-level subsystem become correlated to orthogonal continuum wave packets. 

The dissipative and unitary tunneling rates Rd and R u can also be obtained from the 
reduced density in a straightforward way. A typical result for moderately strong coupling is 
shown in fig. 3. 

3.4 Interpretation of model results. 

In spite of the schematic character of the model, the features emerging from its analysis 
underline relevant aspects which should be kept in mind in the study of more realistic cases. 
As shown explicitly e.g. in the general equation flT7|), the coherence properties directly 
relevant for localization are, in the representation used there, contained in the real part of 
the matrix element p+-. The harmonic time dependence of this quantity in the absence of 
any external entanglements is modified when these effects are turned on by making g and 
g' (or, equivalently, T and T') different from zero, as can be read from the time dependent 
term of equation (|26|) . In the case considered in section [3.2| (T ^ 0, V = 0) the harmonic 
oscillations are damped by the effects of the external coupling, leading, in the strong coupling 
limit, to tunneling rates dominated by the associated relaxation times, rather that by the 
oscillation period. When also T' 7^ 0, on the other hand, a new undamped oscillatory term 

arises due to contributions to Re p H coming from continuum components, an effect which 

has been referred to as a re-correlation process in section The amplitude of the damped 
component is at the same time reduced, leading to its complete cancellation in the limiting 
case r = r" 7^ 0. 

It is important to stress that the coherence properties bearing on the purity of the 
complete tunneling state, which can be measured in a basis-independent way in terms of 
the determinant of the reduced density matrix (see equation (|T8|)), involve other dynamic 

quantities besides the localization-specific correlation part Rep H , and therefore in general 

correlate poorly with the tunneling probability P{t). In particular, when T = V ^ 
the transient decoherence undergone by the tunneling state has no effect on the tunneling 
probability, since the localization-specific coherence properties remain the same as in the 
absence of the external couplings due to the re-correlation process. The decomposition of 
the tunneling rate P(t) into dissipative and unitary parts Rd{t) and R u (t) can thus be seen as 
expressing the dependence of the time rate of change of the localization-specific correlation 
part Re p+_ in terms of contributions associated to the (unitary) time evolution of the 




1 



2rr 



2 
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eigenvectors and from the (non-unitary) time evolution of the eigenvalues of the reduced 
density which describes the tunneling state. Only the latter contributions relate to the 
overall decoherence of the tunneling state. 



4 Master equation results. 

The dissipative dynamics of the reduced density matrix of a two-level system is often de- 
scribed in terms of a master equation valid in a weak-coupling regime, derived under the 
assumption of linear coupling to a reservoir of harmonic oscillators, which is subsequently 
eliminated from the description under the assumption of the validity of a Born-Markov ap- 
proximation. In addition, one assumes in the derivation that correlation functions involving 
reservoir degrees of freedom can be meaningfully evaluated for all times using the initial 
state of the reservoir, usually taken to be a thermal equilibrium state. The usual form of 
the master equation || is obtained assuming that the coupling to the reservoir is given by 

H mt = Y (9kbl<T- + 9* k ha^ 

k 

where b k , b\ are the lowering and rising operators for the fc-th bath oscillator. This coupling 
term is similar to the term involving g(rj) in eq. @ if one associates the ground state of 
the reservoir with the discrete state |0&) of the model. The resulting form of the equation is 
then, in the limit of zero temperature, 

P(t) = -«| fo.pC*)] + | (2a_p(t)o+ - <r+<r-p(t) - p(t)a + a-) 

where 7 is a damping coefficient related to H mt and to the distribution of reservoir frequencies. 
The solution of this equation satisfying the initial condition 

P++ (0)= P+ _(0) = ^ (27) 



is 



2 



.(f) = \e-* and P+ _(t) = 



which reproduces the result (0), with the correspondence 7 T, up to details related 
to the shift function F(E) in that case. It should be kept in mind, however, that the 
derivation of the master equation is limited to weak coupling regimes due to the Born- 
Markov approximation. 

One may next try and include effects analogous to those produced by the coupling term 
involving g'(rj) in the model Hamiltonian ([D by adding a new term to H int in the derivation 
of the master equation. Thus, using now 

Hmt -> Yl (9iA<T- + g* k b k cr + ) + (9k b k°+ + g'lbkV- 
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we obtain, following the standard derivation procedure 0, the modified master equation 



pit) = ~i\ [ ff 3, pit)] + \ i2v-p(t)a + - a + a_p(t) - p(t)a + a_) (28) 

+ 1 - (2a + p(t)a_ - 0-a + p(t) - p(t)a_o + ) + ^ (2a +P (t)a + + 2a_p(t)a_) 

where 7' a the constant analogous to 7 but related to the new terms added to H int . The 
solution of this equation satisfying the same initial condition ( p7| ) is 

p++(f) = i e -w>< + ^(i- e -wv><) 

and 

p-(7+7')* 



77' — ie sen (2\Je 2 — 77' + \Ze 2 — 77' cos (2\Je 2 — 77' 



2^ - 77' 

which is different from the corresponding results obtained in subsection |3.3| . In particular, 
the model result for p ++ (t), eq. (p4|), involves two relaxation times associated to the damp- 
ing widths T and T' appearing in the difference of separate exponentials, while the master 
equation result involves just one "effective" relaxation time appearing in a single exponen- 
tial involving the sum 7 + 7'. The relation between these two results can be revealed by 
expanding both solutions to first order in the damping parameters. One obtains, from the 
master equation result, 

which is identical to what one obtains from eq. (^4j) with the replacements 7 <->• V and 
7' <-» r'. The expansion of the master equation result for p+-(t) coincides likewise with the 
expansion of eq. (|25|). These facts can be understood by recalling that a Born-Markov ap- 
proximation is involved in the derivation of the master equation. The Born approximation 
implies a truncation of an expansion in the couplings, producing a result which is subse- 
quently "exponentiated" to all orders through the Markov assumption. As a result of this 
there is no general guarantee for results of the master equation which go beyond the order 
of the truncation involved in the Born approximation. In particular, the dynamics leading 
to two different relaxation times is not properly accounted for. Note, in this connection, 
that when 7 = 7' (i.e. when the two relaxation times coincide) the master equation result 
reproduces the model result up to a small shift of the bare frequency e, if one restricts oneself 
to weak coupling regimes. 

It may be argued that the master equation (p§|) is derived on the assumption of a reservoir 
of harmonic oscillators, which differs from the dynamical assumptions of the model defined by 
the Hamiltonian (§) . However, the same master equation also follows from this Hamiltonian, 
with the usual (in fact, even somewhat weaker) derivation assumptions, establishing, in 
particular, the correspondence 7 <-> T and 7' <-> V. The Born approximation has to be 
performed in the same way in both cases, indicating again that it is at the root of the 
discussed discrepancies. 
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5 Concluding remarks. 



The main general conclusion to be drawn from the analysis of the model described by the 
Hamiltonian (|S|) with localized initial condition ([TJ]) is that a clear distinction should be 
made between coherence properties relating a) to localization in the tunneling degree of 
freedom and b) to the degree of quantum purity of the reduced density (§) describing the 
quantum state of the tunneling system, as measured e.g. by the linear entropy (||). Coherence 
properties of class a) are carried by the real part of off-diagonal matrix element p+_ (t) in the 
basis which diagonalizes H a or, equivalently, by the diagonal matrix elements in the localized 
basis {|r), e.g. p rr (t) = 1/2 — Rep + _ = P(t). Coherence properties of class b), on the 
other hand, are carried by the eigenvalues of the reduced density, and involve therefore other 
dynamic quantities in addition to that which is relevant for class a). 

The contribution of decoherence processes undergone by the reduced density to the coher- 
ence properties of class a) can be isolated through the general decomposition of the tunneling 
rate P(t) into dissipative and unitary parts, eqs. (0) and (H). Results for the damped two- 
level model show that the dissipative component contributes most at times short on the scale 
of the bare frequency e. The enhanced initial tunneling rate in the strong coupling case with 
no re-correlation effects (see fig. 0) is dominated by the dissipative component, and can thus 
be interpreted in terms of a large effect of the coherence loss of the tunneling state also on 
the localization-specific correlation properties of the initial state of the two-level subsystem. 

Comparison of the exact solution of the model with master equation results derived 
in the usual way from the model Hamiltonian indicates that master equations appear to 
be inadequate to deal with situations involving more than a single relaxation time. This 
difficulty can be ascribed to the Born approximation involved in the derivation of these 
equations. 
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Appendix A. 

A general reduced density matrix can always be written in terms of its time dependent 
natural orbitals \n(t)) and associated eigenvalues p n {t) as 

p(t)=J2Ht))p n (t)(n(t)\. 

n 

The probability that the state thus described is found in a given (time independent) subspace 
defined in terms of a projector 

k 
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where the \bk) constitute an orthonormal set of state vectors, is given by 

P(t) = Tr [p(t)V] . 

The vectors \bk) can be expanded in the natural orbitals of pit) as \bk) = J2n PnKt)\ n (t)) so 
that P(t) is expressed as 

P(t)=£pn(*) £|/ff. 
n k 

Derivation with respect to time then gives dissipative and unitary contributions proportional 
to p n and (3^ respectively. 



Appendix B. 

In view of the widespread use of coordinate coupling to a collection of harmonic oscillators in 
order to introduce damping effects on the dynamics of simple quantum-mechanical systems 
(including two- level systems), and since the results obtained in this way seem to be at 
variance with those of the continuum model for the damping as used in section [3], this 
Appendix deals briefly with a comparison of features of the two models. Instead of the 
implementation which minimizes "threshold" effects adopted in that section, the case eo — 
rjo ~ with g(r]) = g*(i]) = g'(f]) of the continuum model will be considered, together with 
with the "zero bias" spin-boson model [[|. 

In both cases the Hamiltonian can be split as H a + H b + H int with H a = |o"3 and 

H b = \0 b )e (0 b \ + [ V dr)\r))r](r)\, H int = a l f drj g(r]) (\rj)(O b \ + \0 b )(r)\) = (?iG 

for the continuum model of equation (S) and 



H b = Y^ Twu a a) a a a , H int = a% £ g a J — (a* + a a ) = o x G. 

for the spin-boson model. A useful way of representing these operators uses the spin basis of 
localized states in which o\ is diagonal. The Hamiltonian is then represented in both cases 
as 



H 



H b + G — |e 



|e H b — G 



The operators H b ± G can be separately diagonalized exactly in the two cases. This solves 
the problem completely in the extreme adiabatic limit e — >• 0. In general, matrix elements in 
the off-diagonal blocks involving H a are modified by overlap factors involving eigenstates of 
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the two operators H b ± G. Rather than dealing with the full problem, we restrict ourselves 
to the perturbative effects of for the lowest adiabatic states. These are, in both cases, 
a degenerate doublet which we denote by \4>q), with energy E . The secular matrix to be 
considered is then 

' E -i e (0(+)|0(-)> \ 

In the case of the spin-boson model the overlap (0o + ^l0o ^) involves oscillator coherent 
states with opposite displacements which depend, in particular, on the coupling constants g a . 
The ensuing reduction of the overlap implies the reduction of the splitting of the perturbed 
states, favoring the slowing down of transition rates out of localized states. In the case of 
the continuum model, on the other hand, E is the root which lies below r] of the dispersion 
equation 

E -e = dr] — . 

J m E -r) 

The corresponding unperturbed states can be expanded as 

\4 ±) ) = <& ±) \0 b ) + fdri A^\n)\n) 
the expansion coefficients being given by 

|a^| 2 = — * ^- and A^\rj) = . 



and the overlap factor is 



_ i rv j n 9 2 (v) 

i.(+)u(-)\ _ 4 J V0 a 'l (Ep-r,) 2 

(v) 



-v) 2 



The behavior of this object involves the energy dependence of the coupling function g{rj) and 
of the solution E of the dispersion equation. In the case when g 2 (rf) = 777 s , with < s < 1 
one finds, for fixed eo = t] = that E approaches minus infinity and |ao ± ' ) | 2 approaches 
unity as the high energy cut-off fj increases. In the same limit the overlap also approaches 
unity and the transition rate out of the localized states approaches its free value. For values 
of fj not too large in comparison with e the overlap is smaller than one, leading to a quenching 
of this transition rate. These behaviors are supported by exact numerical results. 
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Figure captions. 



Fig. 1 - Tunneling probability P(t), idempotency defect S(t) (top) and tunneling rates Rd(t), 
R u (t) and P(t) = Rd(t) + R u (t) (bottom) in the single finite relaxation time, weak-coupling 
case with parameters T = 0.3 and e = 1. The time scale is defined so that h — 1. 

Fig. 2 - Same as Figure 1 in the single finite relaxation time, strong coupling case with 
parameters T = 3., e = 1. Note the smaller depth of the time scale in this case. 

Fig. 3 - Same as Figure 1 in the case of two finite relaxation times with parameters T = 2, 
r = .5, e = 1. 
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Figure 1: Tunneling probability P(t), 
idempotency defect 8(t) (top) and tunnel- 
ing rates Rd(t), R u (t) and P(t) = Rd(t) + 
Ru{t) (bottom) in the single finite relax- 
ation time, weak-coupling case with pa- 
rameters T = 0.3 and e = 1. The time 
scale is defined so that U — 1. 
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Figure 2: Same as Figure |l] in the of sin- 
gle finite relaxation time, strong coupling 
case with parameters T = 3., e = 1. Note 
the smaller depth of the time scale in this 
case. 
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Figure 3: Same as Figure [I] in the case of 
two finite relaxation times, with parame- 
ters r = 2, V = .5, e = 1. 
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